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Introduction 


The  present  study  is  based  on  significant  progress  in  density  functional  theory  (DFT),  and  associated 
software  technology,  which  is  sufficiently  mature  for  the  determination  of  dielectric  response  structure,  and 
should  actually  provide  complementary  information  to  that  obtained  from  experiment.  Calculations  are 
presented  of  excited  state  resonance  structure  associated  with  Fe,  Mn  and  Mg  water  complexes  using  time- 
dependent  density  functional  theory  (TD-DFT),  which  is  an  extension  of  DFT.  Calculation  of  excited  state 
resonance  structure  using  TD-DFT  can  provide  interpretation  of  absorption  spectra  with  respect  to  molecular 
structure  for  excitation  by  electromagnetic  waves  at  frequencies  within  the  UV  range.  In  principle,  these 
absorption  spectra  should  provide  quantitative  initial  estimates  of  spectral  response  features  that  can  be 
subsequently  adjusted  with  respect  to  additional  information  such  as  laboratory  measurements  and  other 
types  of  theory  based  calculations,  or  conversely,  adapted  as  constraints  for  the  inverse  analysis  of 
experimentally  measured  absorption  spectra.  A  significant  aspect  of  using  DFT  and  TD-DFT  for  the 
calculation  of  absorption  spectra  is  that  it  adopts  the  perspective  of  computational  physics,  according  to 
which  a  numerical  simulation  represents  another  source  of  “experimental”  data. 

Density  functional  theory  has  been  successfully  used  to  investigate  the  vibrational  spectra  of 
energetic  materials  in  the  form  of  single  molecules  and  molecular  crystals  [1-9].  These  calculations  provide 
interpretation  of  absorption  spectra  with  respect  to  molecular  structure  for  various  forms  of  materials,  which 
can  be  encountered  in  various  detection  scenarios.  In  particular,  the  calculated  absorption  spectra  of  isolated 
molecules  can  help  to  identify  intramolecular  vibrational  modes  of  various  materials.  A  series  of  studies  have 
focused  on  the  general  concept  of  constructing  dielectric  response  functions  using  DFT  for  the  purpose  of 
quantitative  simulation  of  explosives  detection  scenarios  [9,10,11].  For  these  studies  the  DFT  software 
GAUSSIAN09  (G09)  was  adopted  [12].  As  emphasized  in  these  studies,  the  calculation  of  absorption  spectra 
using  DFT  and  TD-DFTdefines  a  general  approach  where  dielectric  response  is  estimated  within  the  bounds 
of  relatively  well-defined  adjustable  parameters.  Following  this  approach,  permittivity  functions  are 
constructed  using  DFT  and  TD-DFT  calculated  absorption  spectra  under  the  condition  that  the  calculated 
resonance  locations  are  fixed,  while  resonance  widths  and  number  densities  are  assumed  adjustable  with 
respect  to  additional  information  such  as  experimentally  observed  spectra  or  more  advanced  theory. 

Previous  studies  have  examined  various  properties  of  water  molecules  and  their  clusters  [13-21],  The 
absorption  spectrum  of  molecular  clusters  consisting  of  water  molecules  should  be  of  significance  for 
interpretation  of  absorption  spectra  associated  with  detection  in  practice.  This  follows  in  that  most 
environments  associated  with  detection  in  practice  include  the  presence  of  water  in  one  form  or  another. 
These  forms  can  range  from  isolated  molecules  in  gas  phase,  molecular  clusters,  adsorbed  surface  layers, 
droplets  and  interface  regions  in  liquid  phase,  and  ice.  Absorption  spectra  of  metal-water  complexes 
represent  an  important  case  for  dielectric  response  with  respect  to  electromagnetic  wave  excitation.  This 
regime  should  be  better  quantified  for  improved  interpretation  of  absorption  spectra  associated  with  systems 
that  include  metal-water  complexes  as  components. 

The  organization  of  the  subject  areas  presented  are  as  follows.  First,  a  general  review  of  the  elements 
of  vibrational  analysis  using  DFT  for  determination  the  stability  of  molecular  geometies  is  presented. 
Second,  a  brief  discussion  presedented  concerning  the  significance  of  imaginary  frequencies  for  vibrational 
analysis  using  DFT.  Third,  the  excited  state  resonance  structure  of  Fe,  Mn  and  Mg  water  complexes  is 
investigated  using  TD-DFT  calculations.  Finally,  a  conclusion  is  given. 


Manuscript  approved  July  1 1,  2013. 


i 


Calculation  of  Absorption  Spectra  using  DFT 


As  in  previous  studies  [9,  10,  11]  the  formal  mathematical  structure  underlying  DFT  calculations,  as 
well  as  the  procedure  for  calculation  of  absorption  spectra  corresponding  to  vibrational  states,  is  included 
here  for  purposes  of  completeness.  The  extension  of  DFT  for  the  calculation  of  absorption  spectra 
corresponding  to  electronic  excitation  states,  which  is  the  formalism  of  time-dependent  density  functional 
theory  (TD-DFT),  is  described  in  reference  [22],  In  this  study,  we  consider  the  calculation  of  absorption 
spectra  for  electronic  excitation  states.  The  calculation  of  absorption  spectra  for  vibrational  states  is 
considered,  however,  for  the  purpose  of  investigating  the  stability  of  molecular  geometries.  With  respect  to 
the  relative  stability  of  molecular  geometries,  the  significance  of  imaginary  frequencies  is  discussed  below. 

The  DFT  software  GAUSSIAN09  (G09)  can  be  used  to  compute  an  approximation  of  the  IR 
absorption  spectrum  of  a  molecule  or  molecules  [23].  This  program  calculates  vibrational  frequencies  by 
determining  second  derivatives  of  the  energy  with  respect  to  the  Cartesian  nuclear  coordinates,  and  then 
transforming  to  mass-weighted  coordinates  at  a  stationary  point  of  the  geometry.  The  IR  absorption  spectrum 
is  obtained  using  density  functional  theory  to  compute  the  ground  state  electronic  structure  in  the  Bom- 
Oppenheimer  approximation  using  Kohn-Sham  density  functional  theory  [24-29].  GAUSSIAN  uses 
specified  orbital  basis  functions  to  describe  the  electronic  wavefunctions  and  density.  For  a  given  set  of 
nuclear  positions,  the  calculation  directly  gives  the  electronic  charge  density  of  the  molecule,  the  potential 
energy  V,  and  the  displacements  in  Cartesian  coordinates  of  each  atom.  The  procedure  for  vibrational 
analysis  followed  in  GAUSSIAN  is  that  described  in  Ref  [23].  Reference  [29]  presents  a  fairly  detailed 
review  of  this  procedure.  A  brief  description  of  this  procedure  is  as  follows. 

The  procedure  followed  by  GAUSSIAN  is  based  on  the  fact  the  vibrational  spectrum  depends  on  the 
Hessian  matrix  fcART,  which  is  constructed  using  the  second  partial  derivatives  of  the  potential  energy  V  with 
respect  to  displacements  of  the  atoms  in  Cartesian  coordinates.  Accordingly,  the  elements  of  the  3N  x  3N 
matrix:  fcART  are  given  by 
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coordinates,  and  N  is  the  number  of  atoms.  As  discussed  above,  the  zero  subscript  in  Eq.(l)  indicates  that  the 
derivatives  are  taken  at  the  equilibrium  positions  of  the  atoms,  and  that  the  first  derivatives  are  zero.  Given 
the  Hessian  matrix  defined  by  Eq.(l)  the  operations  for  calculation  of  the  vibrational  spectrum  require  that 
the  Hessian  matrix  Eq.(l)  be  transformed  to  mass-weighted  Cartesian  coordinates  according  to  the  relation 
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weighted  Cartesian  coordianates.  GAUSSIAN  computes  the  energy  second  derivatives  Eq.(2),  thus 
computing  the  forces  for  displacement  perturbations  of  each  atom  along  each  Cartesian  direction.  The  first 
derivatives  of  the  dipole  moment  with  respect  to  atomic  positions  d/j,  /<?§,.  are  also  computed.  Each  vibrational 
eigenmode  leads  to  one  peak  in  the  absorption  spectrum,  at  a  frequency  equal  to  the  mode's 
eigenfirequencyv/;0.  The  absorption  intensity  corresponding  to  a  particular  eigenmode  n  whose 
eigenfrequency  is  vw0  is  given  by 
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where  Icart  is  the  matrix  whose  elements  are  the  displacements  of  the  atoms  in  Cartesian  coordinates.  The 
matrix  Icart  is  determined  by  the  following  procedure.  First, 


ICART  -  MImwC  , 


(Eq4) 


where  Imwc  is  the  matrix  whose  elements  are  the  displacements  of  the  atoms  in  mass-weighted  Cartesian 
coordinates  and  M  is  a  diagonal  matrix  defined  by  the  elements 
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(Eq5) 


Proceeding,  lMwc  is  the  matrix  needed  to  diagonalize  fMwc  defined  by  Eq.(2)  such  that 

(  1mwc)T  fMwc(  Imwc)  =  A  ,  (Eq  6) 

where  A  is  the  diagonal  matrix  with  eigenvalues  A,,  The  procedure  for  diagonalizing  Eq.(6)  consists  of  the 
operations 

fiNT  =  (  D)T  fMwc(  D)  (Eq  7) 

and 

(  L)t  fMwc(  L)  =  A  ,  (Eq  8) 


where  D  is  a  matrix  transformation  to  coordinates  where  rotation  and  translation  have  been  separated  out  and 
L  is  the  transformation  matrix  composed  of  eigenvectors  calculated  according  to  Eq.(8).  The 
eigenfrequencies  in  units  of  (cm1)  are  calculated  using  the  eigenvalues  kn  by  the  expression 
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where  c  is  the  speed  of  light.  The  elements  of  Icart  are  given  by 
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where  k,  i=l,...,  3N,  and  the  column  vectors  of  these  elements  are  the  normal  modes  in  Cartesian 
coordinates. 

The  intensity  Eq.(3)  must  then  be  multiplied  by  the  number  density  of  molecules  to  give  an 
absorption-line  intensity  in  the  non-interacting  molecule  approximation.  It  follows  that  the  absorption 
spectrum  calculated  by  GAUSSIAN  is  a  sum  of  delta  functions,  whose  line  positions  and  coefficients 
correspond  to  the  vibrational-transition  frequencies  and  the  absorption-line  intensities,  respectively.  In 
principle,  however,  these  spectral  components  must  be  broadened  and  shifted  to  account  for  anharmonic 
effects  such  as  finite  mode  lifetimes  and  inter-mode  couplings. 

The  DFT  software  GAUSSIAN09  (G09)  can  be  used  to  compute  an  approximation  of  the  UV 
absorption  spectrum  of  a  molecule  or  molecules  by  application  of  the  formalism  of  time-dependent  density 
functional  theory  (TD-DFT),  which  is  described  in  reference  [22], 
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Physical  Significance  of  Quantum  Mechanical  Transition  State  for 
Molecular  Transformation 


A  molecule  in  3-dimensions  has  a  total  of  3N-6  normal  mode  vibrations.  The  Schrodinger  equation  for 
the  harmonic  oscillations  of  these  normal  modes  has  known  solutions.  The  quantum  mechanical  spectrum  of 
each  of  these  vibrations  is  given  in  the  harmonic  approximation  by  the  energies  En  =  (n+l/2)h  v,  where  n  is  a 
quantum  number,  h  is  Planck’s  constant,  and  v  is  a  vibration  frequency  given  by  v=V(k/m),  where  k  is  the 
spring  constant  of  the  normal  vibration  and  m  is  the  effective  mass  contributed  by  those  atoms  vibrating  in 
the  normal  mode.  A  molecule  in  stable  equilibrium  is  characterized  by  all  positive  normal  mode 
frequencies  v.  But  the  definitive  mathematical  characteristic  of  a  transition  state  is  that  it  has  all  positive 
frequencies  but  one,  which  is  imaginary.  That  is  to  say,  (k/m)  is  a  negative  number.  The  vibration 
corresponding  to  an  imaginary  frequency  is  one  in  which  the  atoms  are  breaking  away  from  bonds 
characteristic  of  chemical  reactants  and  are  moving  towards  those  bonds  characteristic  of  chemical  products. 
Chemical  reactions  break  bonds  in  the  reactants,  rearrange  them  and  form  new  bonds  in  the  products.  The 
transition  state  is  a  particular  geometric  arrangement  of  the  atoms  in  a  chemical  system,  at  the  maximal  peak 
of  the  energy  surface  separating  reactants  from  products.  In  the  transition  state  every  normal  vibration 
distortion  but  one,  occurs  within  a  stable  energy  minimum.  But  the  one  normal  mode  distortion  of  imaginary 
frequency  occurs  at  an  unstable  energy  maximum  sending  reactants  toward  products.  The  height  of  energy 
peak,  the  activation  energy  Ea,  associated  with  the  occurrence  of  a  transition  state  determines  the  minimal 
energy  accumulated  by  reactants  to  surmount  the  barrier  separating  reactants  from  products.  In  the  Arrhenius 
formulation,  the  rate  constant  for  the  reaction  is  given  by  T  =  A  e  “Ea/RT,  where  A  is  a  constant,  R  is  the  ideal 
gas  constant,  and  T  is  the  absolute  temperature.  The  common  occurrence  is  that  a  particular  transition  state 
mechanism  for  chemical  reaction  is  associated  with  one  imaginary  frequency,  and  therefore  a  single 
mechanism  of  reaction.  Much  less  commonly,  an  energy  surface  of  multiple  channels  of  reaction  mechanism 
may  give  rise  to  correspondingly  multiple  imaginary  frequencies.  But  in  any  event,  the  transition  state 
contains  the  energetic  and  geometric  information  that  defines  the  transformation  inherent  within  chemical 
reactions. 

Excited  State  Resonance  Structure  of  Fe,  Mn  and  Mg  Water  Complexes 

In  this  section  are  presented  the  results  of  computational  investigations  using  DFT  and  TD-DFT 
concerning  Fe,  Mn  and  Mg  water  complexes  of  various  size.  These  results  include  the  relaxed  or  equilibrium 
configurations  of  metal-water  complexes,  denumeration  of  imaginary  frequencies  for  stability  analysis  of 
molecular  geometries,  which  are  calculated  by  DFT,  and  oscillator  strength  as  a  function  of  excitation  energy 
(within  the  UV  range)  for  different  geometries  of  the  interacting  systems  associated  with  stable  structures, 
which  are  calculated  by  TD-DFT.  For  these  calculations  geometry  optimization  and  vibrational  analysis  was 
effected  using  the  DFT  model  B3LYP  [30,  31]  and  basis  function  6-31G(d,p)  [32,  33].  According  to  the 
specification  of  this  basis  function,  (d,p)  designates  polarization  functions  having  one  set  of  d  and  p 
functions  for  heavy  atoms  and  hydrogen  atoms,  respectively  [34], 

A  graphical  representation  of  Fe,  Mn,  Mg  water  complexes  having  different  molecular  geometries  and 
sizes  are  shown  in  Figs.  (1)  through  (25).  Given  in  Tables  1  through  5  are  excitation  energies  and  oscillator 
strenghts  for  various  molecular  geometries,  sizes,  and  electronic  structures  of  Fe,  Mn,  Mg  water  complexes 
with  and  without  water  background  after  geometry  optimization  calculated  using  TD-DFT.  Given  in  Table  6 
are  total  energies  for  optimized  geometries  and  denumeration  of  imaginary  frequencies  for  calculation  of 
vibrational  (IR)  spectra  using  DFT  (B3LYP/6-31G(d,p))  for  Mn  water  complexes.  Shown  in  the  appendix 
are  the  geometrie  of  Fe++-6H20  (D3d)  and  (D4h)  corresponding  to  stable  structures  (see  Figs.  (Al)  and  (A2)). 
For  these  geometries  TD-DFT  calculations  predict  that  there  are  no  UV  excitations. 
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Conclusion 


The  TD-DFT  calculated  absorption  spectra  given  here  provide  two  types  of  information  for  general 
analysis  of  dielectric  response.  These  are  the  denumeration  of  excited  state  resonance  modes  and  estimates  of 
molecular  level  dielectric  response  structure.  The  calculations  of  excited  state  resonance  structure  associated 
with  Fe,  Mn  and  Mg  water  complexes  using  TD-DFT  are  meant  to  serve  as  reasonable  estimates  of 
molecular  level  response  characteristics,  providing  interpretation  of  dielectric  response  features,  for 
subsequent  adjustment  relative  to  experimental  measurements  and  additional  constraints  based  on  molecular 
structure  theory.  With  respect  to  spectroscopic  methods  for  detection  or  monitoring  of  target  molecules 
within  various  types  of  water  environments,  i.e.,  different  types  of  detection  strategies  and  their  associated 
algorithms  for  post-processing  of  measurements,  the  calculated  resonance  spectra  presented  here  serve  the 
purpose  of  isolating  spectral  features  for  correlating  with  the  presence  of  metal  water  complexes. 
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Appendix.  Stable  geometries  of  Fe++  -6(H20)  showing  no  UV  excitations  according  to  TD-DFT 


Figure  Al.  Molecular  geometry  of  Fe++-6H20  (D3<j)  corresponding  to  stable  structure.  For  this  geometry 
there  is  no  UV  excitation. 


Figure  A2.  Molecular  geometry  of  Fe++-6H20  (D4h)  corresponding  to  stable  structure.  For  this  geometry 
there  is  no  UV  excitation. 
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Table  1.  Excitation  energies  and  oscillator  strenghts  for  Fe++  -3(H20)  and  Fe  *3(H20)  after 


geomel 

;ry  optimization  calculated  using  TD-DFT. 

Complex 

Excitation 

Energy 

Oscillator 

Strength 

Complex 

Excitation 

Energy 

Oscillator 

Strength 

Fe++  *3(H20) 

847.83  nm 
831.67  nm 
667.04  nm 

0.0003 

0.0004 

0.0005 

Fe++  *3(H20) 
in  water 

1958.26  nm 
1560.28  nm 
762.19  nm 
707.94  nm 
641.34  nm 

0.0001 

0.0001 

0.0002 

0.0002 

0.0003 

Fe  *3(H20) 

609.82  nm 

0.0326 

Fe  *3(H20) 
in  water 

637.73  nm 

0.0356 

Table  2.  Excitation  energies  and  oscillator  strenghts  for  Fe++  -5(H20)  after  geometry  optimization  calculated 


using  TD-DFT  (] 

B3LYP/6-31G(d,] 

P))- 

Complex 

Excitation 

Energy 

Oscillator 

Strength 

Complex 

Excitation 

Energy 

Oscillator 

Strength 

Fe++  5(H20) 

0>3d) 

2007.74  nm 
514.92  nm 
484.82  nm 

0.0001 

0.0004 

0.0004 

Fe+2  *5(H20) 

0>3d) 

In  water 

1473.84  nm 
1322.75  nm 
489.22  nm 
467.55  nm 

0.0002 

0.0001 

0.0002 

0.0003 

Figure  1.  Molecular  geometry  of  Fe++-3H20  corresponding  to  stable  structure. 


Figure  2.  Molecular  geometry  of  Fe-3H20  corresponding  to  stable  structure. 
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Figure  3.  Molecular  geometry  of  Fe++  -5(H20)  (D3d)  corresponding  to  stable  structure. 


_i_2 

Figure  4.  Molecular  geometry  of  Mn  •  H20  corresponding  to  stable  structure.  For  both  doublet  and  singlet 
multiplicity  there  is  no  UV  excitation. 


1 2 

Figure  5.  Molecular  geometry  of  Mn  •  3H20  corresponding  to  stable  structure.  For  the  case  of  doublet 
multiplicity  there  is  no  UV  excitation. 
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I  o 

Figure  6.  Molecular  geometry  of  Mn  •  4H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 


|  O 

Figure  7.  Molecular  geometry  of  Mn  •  4H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 
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Figure  8.  Molecular  geometry  of  Mn+3  •  6H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 
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I  o 

Figure  9.  Molecular  geometry  of  Mn  •  6H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 


I  O 

Figure  10.  Molecular  geometry  of  Mn  •  6H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  no  UV  excitation. 
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Figure  11.  Molecular  geometry  of  Mn+3  •  7H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 
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I  o 

Figure  12.  Molecular  geometry  of  Mn  •  7H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 


Figure  13.  Molecular  geometry  of  Mn  •  8H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  UV  excitation. 


Figure  14.  Molecular  geometry  of  Mn+3  •  8H20  corresponding  to  stable  structure.  For  the  case  of  singlet 
multiplicity  there  is  no  UV  excitation. 
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_i_2  1 2 

Table  3.  Excitation  energies  and  oscillator  strenghts  for  Mn  •  n(H20)  (doublet)  and  Mn  •  n(H20) 
(singlet)  after  geometry  optimization  calculated  using  TD-DFT  (B3LYP/6-31G(d,p)). 


Mn+2  •  n(H20)  (doublet) 

Mn+3  •  n(H20)  (singlet) 

n(H20) 

Excitation  E 

Oscillator  strength 

n(H20) 

Excitation  E 

Oscillator  strength 

n  =1 

no  peak 

Oscillator  strength  =  0 

n  =  2 

Unstable  geometry 

Unstable  geometry 

n  =  3 

no  peak 

Oscillator  strength  =  0 

n  =  4 

None 

(Convergence  criterion  not 
met) 

n  =  5 

None 

(Convergence  criterion  not 
met) 

n  =  6 

Unstable  geometry 

n  =  7 

None 

(Convergence  criterion  not 
met) 

n  =  8 

None 

(Convergence  criterion  not 
met) 

n  =  12 

None 

(Convergence  criterion  not 
met) 

n  =  1 

no  peak 

Oscillator  strength  =  0 

n  =  2 

Unstable  geometry 

n  =  3 

Unstable  geometry 

348.11  nm 

0.0137 

11 12.09  nm 

0.0002 

n  =  4 

1111.13  nm 

0.0002 

773.28  nm 

0.0016 

625.73  nm 

0.0024 

625.08  nm 

0.0024 

Unstable  geometry 

Unstable  geometry 

Unstable  geometry 

Unstable  geometry 

762.69  nm 

0.0007 

759.70  nm 

0.0012 

n  —  o 

717.17  nm 

0.0022 

514.02  nm 

0.0029 

no  peak 

Oscillator  strength  =  0 

no  peak 

Oscillator  strength  =  0 

904.18  nm 

0.0001 

751.21  nm 

0.0023 

671.68  nm 

0.0014 

n  =  7 

542.40  nm 

0.0036 

1 177.44  nm 

0.0002 

982.24  nm 

0.0017 

614.12  nm 

0.0013 

561.69  nm 

0.0013 

no  peak 

Oscillator  strength  =  0 

n  =  8 

Unstable  geometry 

n  =  12 

Unstable  geometry 
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Table  4.  Excitation  energies  and  oscillator  strenghts  for  Mg  •  n(H20)  with  and  without  water  background 
after  geometry  optimization  calculated  using  TD-DFT  (B3LYP/6-31G(d,p)). 


Componds 

Excitation  E 

Oscillator 

strength 

Componds 

Excitation  E 

Oscillator 

strength 

Mg+2  1(H20) 

233.27  nm 
136.24  nm 
124.18  nm 
120.26  nm 
114.83  nm 

0.0080 

0.0590 

0.1341 

0.0409 

0.0110 

Mg+2  1(H20)  in 
water 

157.32  nm 
124.55  nm 

1 12.92  nm 
109.63  nm 

0.0288 

0.0156 

0.0020 

0.1892 

Mg+2  -2(H20) 

175.98  nm 
175.98  nm 
124.41  nm 

0.0201 

0.0201 

0.0079 

Mg+2  -2(1 120)  in 
water 

152.17  nm 

152.17  nm 

127.18  nm 
127.18  nm 

0.0341 

0.0341 

0.0117 

0.0117 

Mg+2  3(H20) 

152.88  nm 
152.14  nm 
121.01  nm 

0.0303 

0.0588 

0.0432 

Mg+2  -3(H20)  in 
water 

147.99  nm 
147.35  nm 
127.31  nm 
126.89  nm 

0.0443 

0.0732 

0.0189 

0.0083 

Mg+2  -4(H20) 

- 

- 

Mg+2  -4(H20)  in 
water 

- 

- 

Mg+2  5(H20) 

146.59  nm 

144.31  nm 

143.31  nm 
142.12  nm 

0.0236 

0.0782 

0.0609 

0.0501 

Mg+2  -5(H20)  in 
water 

148.60  nm 
146.13  nm 
144.94  nm 
144.94  nm 

0.0553 

0.0605 

0.0560 

0.0705 

Mg+2  -6(H20) 

725.79  nm 
144.97  nm 
144.97  nm 
144.97  nm 

0.3477 

0.0940 

0.0940 

0.0940 

Mg+2  -6(H20)  in 
water 

146.83  nm 
146.83  nm 
146.83  nm 

0.1048 

0.1048 

0.1048 
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Table  5.  Excitation  energies  and  oscillator  strenghts  for  Mg-  n(H20)  with  and  without  water 
background  after  geometry  optimization  calculated  using  TD-DFT  (B3LYP/6-31G(d,p)). 


Oscillator 

Oscillator 

Componds 

Excitation  E 

strength 

Componds 

Excitation  E 

strength 

341.28  nm 

0.5217 

373.21  nm 

0.5508 

333.85  nm 

0.5312 

371.06  nm 

0.0029 

Mg  •  1(H20) 

320.04  nm 

0.5599 

Mg  •  1(H20)  in 

362.03  nm 

0.5677 

207.67  nm 

0.1132 

water 

208.95  nm 

0.0870 

198.23  nm 

0.0029 

192.17  nm 

0.0031 

166.73  nm 

0.0299 

172.23  nm 

0.0318 

378.26  nm 

0.6040 

423.86  nm 

0.6413 

372.02  nm 

0.4846 

422.77  nm 

0.5068 

Mg  -2(H20) 

372.02  nm 

0.4846 

Mg  -2(H20)  in 
water 

422.77  nm 

0.5068 

201.22  nm 

0.0079 

198.18  nm 

0.0116 

201.22  nm 

0.0079 

198.18  nm 

0.0116 

Mg  -3(H20) 

- 

- 

Mg  •3(H20)  in 
water 

- 

- 

354.06  nm 

0.5032 

390.51  nm 

0.5549 

335.15  nm 

0.4924 

371.33  nm 

0.4807 

Mg  *4(H20) 

313.65  nm 

0.4655 

Mg  -4(H20)  in 

346.89  nm 

0.4895 

226.00  nm 

0.0377 

water 

234.66  nm 

0.0343 

209.57  nm 

0.0015 

216.34  nm 

0.0089 

198.68  nm 

0.0283 

203.49  nm 

0.0302 

593.22  nm 

0.2000 

517.52  nm 

0.4093 

541.31  nm 

0.4545 

506.00  nm 

0.5147 

Mg  -5(H20) 

535.79  nm 

0.3704 

Mg  -5(H20)  in 

475.78  nm 

0.1890 

327.91  nm 

0.0631 

water 

292.27  nm 

0.0279 

281.44  nm 

0.0693 

259.25  nm 

0.0984 

264.19  nm 

0.0084 

234.97  nm 

0.0005 

725.79  nm 

0.3477 

820.34  nm 

0.3865 

Mg  •6(H20) 

725.79  nm 
725.79  nm 

0.3477 

0.3477 

Mg  •6(H20)  in 
water 

820.34  nm 
820.34  nm 

0.3865 

0.3865 

307.15  nm 

0.0045 

303.45  nm 

0.0037 
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Table  6.  Total  energy  for  optimized  geometry  and  denumeration  of  imaginary  frequencies  for  calculation  of 
vibrational  (IR)  spectra  using  DFT  (B3LYP/6-31G(d,p)). 


Mn+2  •  n(H20)  (doublet) 

Mn+3  •  n(H20)  (singlet) 

Energy opt  (a.u.) 

Imag.  Freq. 

Energy opt  (a.u.) 

Imag.  Freq. 

n  =1 

-1226.400 

0 

n  =  1 

-1225.472 

0 

n  =  2 

-1302.970 

1 

n  =  2 

-1302.139 

1 

-1302.957 

1 

n  =  3 

-1379.465 

0 

n  =  3 

-1378.696 

1 

n  =  4 

None 

(Convergence  criterion  not 
met) 

n  =  4 

-1455.351 

0 

-1455.353 

0 

n  =  5 

None 

(Convergence  criterion  not 
met) 

n  =  5 

-1531.868 

1 

-1531.850 

3 

-1608.374 

2 

-1608.373 

4 

n  =  6 

-1608.982 

1 

n  =  6 

-1608.401 

0 
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Figure  15.  Molecular  geometry  of  Mn 
multiplicity  there  is  no  UV  excitation. 
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I2H2O  corresponding  to  stable  structure.  For  the  case  of  singlet 
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Figure  16.  Molecular  geometry  of  Mg  •  1(H20)  corresponding  to  stable  structure.  For  this  structure  there  is 
UV  excitation. 


Figure  17.  Molecular  geometry  of  Mg+2  ^(F^O)  corresponding  to  stable  structure.  For  this  structure  there  is 
UV  excitation. 
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Figure  18.  Molecular  geometry  of  Mg  -3(1120)  corresponding  to  stable  structure.  For  this  structure  there  is 
UY  excitation. 
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Figure  19.  Molecular  geometry  of  Mg  -5(1120)  corresponding  to  stable  structure.  For  this  structure  there  is 
UY  excitation. 


Figure  20.  Molecular  geometry  of  Mg+2  -6(1120)  corresponding  to  stable  structure.  For  this  structure  there  is 
UY  excitation. 


19 


Figure  21.  Molecular  geometry  of  Mg*  1(H20)  corresponding  to  stable  structure.  For  this  structure 
there  is  UV  excitation. 
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Figure  22.  Molecular  geometry  of  Mg*2(H20)  corresponding  to  stable  structure.  For  this  structure 
there  is  UV  excitation. 
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Figure  23.  Molecular  geometry  of  Mg-4(H20)  corresponding  to  stable  structure.  For  this  structure 
there  is  UV  excitation. 
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Figure  24.  Molecular  geometry  of  Mg-5(H20)  corresponding  to  stable  structure.  For  this  structure 
there  is  UV  excitation. 


Figure  25.  Molecular  geometry  of  Mg-6(H20)  corresponding  to  stable  structure.  For  this  structure 
there  is  UV  excitation. 
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